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Abstract 



o 
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^H ! We use numerical modeling to study the features of parametric (quasi-Cherenkov) 

cooperative radiation arising when an electron bunch passes through a crystal (nat- 
ural or artificial) under the conditions of dynamical diffraction of electromagnetic 
f^ i waves in the presence of shot noise. It is shown that in both Laue and Bragg diffrac- 

tion cases, parametric radiation consists of two strong pulses: one emitted at small 
i-^ ' angles with respect to the particle velocity direction and the other emitted at large 

angles to it. Under Bragg diffraction conditions, the intensity of parametric radia- 
tion emitted at small angles to the particle velocity direction reaches saturation at 
C^ ' sufficiently smaller number of particles than the intensity of parametric radiation 

c/3 . emitted at large angles. Under Laue diffraction conditions, every pulse contains two 

strong peaks, which are associated with the emission of electromagnetic waves at 
the front and back ends of the bunch. The presence of noise causes a chaotic signal 
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c^ 
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jS^ . in the interval between the two peaks. 

1 Introduction 
> 

o 

\Q • The generation of short pulses of electromagnetic radiation is a primary challenge 

^ ■ of modern physics. They find applications in studying molecular dynamics in biological 

• ! objects and charge transfer in nanoelectronic devices, diagnostics of dense plasma and 

Q \ radar detection of fast moving objects. 

cn I The advances in the generation of short pulses of electromagnetic radiation in infrared, 

visible, and ultraviolet ranges of wavelengths are traditionally associated with the devel- 
opment of quantum electronic devices - lasers. Radiation in lasers is generated via induced 

^ ! emission of photons by bound electrons. 

j^ \ Electrovacuum devices, operating in a cooperative regime [HE], have recently become 

considered as an alternative to short-pulse lasers, whose active medium is formed by elec- 
trons bound in atoms and molecules. These are free electron lasers, cyclotron-resonance 
masers, and Cherenkov radiators, whose active medium is formed by electron bunches 
propagating in complex electrodynamical structures (undulators, corrugated waveguides 
and others). The feature of the cooperative operation regime lies in the fact that the 
radiation power scales as the squared number of particles in the bunch. This allows call- 
ing this regime "superradiance" by analogy with the phenomenon predicted by Dicke in 
quantum electronics [Ij. It should be noted, that the electrons moving in FELs are ini- 
tially homogeneously distributed in the phase space. As a result, bremsstrahlung starts 
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from sponteneous emission. This is true even if the bunch length is much smaller than 
the radiation wave length. In contrast to bremsstrahlung, Cherenkov (quasi-Cherenkov) 
radiation starts from coherent emission when such a short-length bunch is injected into 
a slow-wave structure, i. e. the radiation power is proportional to the squared number 
of particles. The question arises whether this dependence holds when the bunch length 
increases. 

This paper considers coperative radiation emitted by electron bunches when charged 
particles pass through crystals (natural or artificial) under the conditions of dynamical 
diffraction of elctromagnetic waves. Note that a detailed analysis of the features of spon- 
taneous radiation of electrons passing through crystals in both frequency [3] and time [1] 
domains has been carried out before. This radiation, emitted at both large and small 
angles with respect to the direction of electron motion, is called the parametric (quasi- 
Cherenkov) radiation. The problems of amplification of induced parametric radiation 
have also been thoroughly studied in the literature [5], and the threshold current densi- 
ties providing lasing in crystals have been calculated [6]. 

This paper is arranged as follows: In the beginning, a nonlinear theory of interaction of 
relativistic charged particles and the electromagnetic field in crystals is set forth, followed 
by the results of numerical calculations of the parametric radiation pulse. Then, the 
dependence of the radiation intensity on the electron bunch length and the geometrical 
parameters of the system is considered. The appendix outlines the algorithm used in 
the simulation. The feature of the algorithm is that it is based on the particle-in-cell 
method, which enables studying kinetic phenomena. Let us note that in most of the 
existing codes (see, e.g. [TirSJ) used for simulating the interaction of charged particles and 
a synchronous wave, the motion of charged particles is considered within the framework 
of the hydrodynamic approximation. 

2 Nonlinear theory of cooperative radiation 

A theoretical analysis of parametric radiation can be performed only by means of a 
self-consistent solution of a nonlinear set of the Vlasov-Maxwell equations: 

| + ,| + ,,(^ + ,x5)| = 0, (1) 

c at 
VxH=-— + -j,j = g,jvfdp, ^2) 

V ■ ^ = Anp, p = qe fd^p, 
V-H = 0, 

describing the electron motion in the electric E and magnetic H fields. Here f{r,p,t) 
is the particle distribution function over the coordinates r and momenta p, while j{r,t) 
and p(r', t) are the current and charge densities, respectively. Since the crystal is a pe- 
riodic linear medium with frequency dispersion, the Fourier transform of the displace- 
ment current D{r,u) relates to the electric field Eif^u) as D{f,uj) = (l + Xoiy) + 
^2xf;(u;) cos(rr))i?(?^, w), where the summation is made over all reciprocal lattice vec- 

tors. Because the dielectric susceptibilities in natural crystals and in grid photonic crystals 



are inversely proportional to the frequency [9]: xo,t{^) = ^or/^'^y Maxwell's equations 
([2]) can be reduced to the equation of the form: 



V(V ■E)-AE + -j-E + 2_, 2^ cos(ff)E = -^. 



C2 9t2 



C2 dt 



(3) 



We shall further be interested in the case when two strong waves are excited in the 
crystal: the forward wave and the diffracted wave (the so-called two- wave diffraction case). 
The forward wave (its wave vector is denoted by k) is emitted at small angles with respect 
to the particle velocity, while the diffracted one, having the wave vector /c,- = fc -f- r, is 
emitted at large angles to it (Fig. 1). Under the conditions of dynamical diffraction, the 
following relation is fulfilled: A;^ ~ /c^ ~ uP' jc? . 






Figure 1: Schemes of parametric radiation in Laue (a, b) and Bragg (c, d) geometries. 



Let us perform the following simplifications: First, we shall neglect the longitudinal 
(V ■ -D — 7- 0) fields of the bunch, which is appropriate when the value of the Langmuir 
oscillation frequency fift of the bunch is less than the values of f2o,r- In this case, the 
Coulomb forces will not have an appreciable effect on dielectric properties of the system. 
Second, we shall seek for the electric field E using the method of slowly varying ampli- 
tudes. Third, we shall assume that a transversally infinite bunch executes one-dimensional 
motion along the OZ-axis (this is achieved by inducing a strong axial magnetic field in 
the system). 

Under the conditions of two-wave diffraction, the field E can be presented as a sum 



E = eoEo{z,t)e 



i(kr—uit) 



+ erET-{z,t)e 



i{krr—idt) 



(4) 



where the amplitudes of the forward Eq and diffracted E^ waves are slowly varying 
variables. This means that for the distances comparable with the wavelength and the 



times comparable with the oscillation period, the values of Eq and E^- remain practically 
the same. Substituting (jlj) into ([1]) and ([3]) and then averaging them over the length 
I = 2jX = 2J7r/kz, where j is a natural number, we obtain 

| + ..| + 2,.Re(^„e'.--.)^^0. (5) 

CZ ZUC ZUIC C J z—l/2 



c dt 



1 dE^ dE^ i^l „ ml „ 27r /'"+'/2 



(72 /cue ZWC C Jz-l/2 



(6) 



Now let us complete the set of equations ([5]) and ([6]) with boundary conditions (the 
initial conditions are reduced to the condition that all values of the fields equal zero at 
i = 0): in the plane z = 0, let us specify the time dependence of function / and set the 
field Eq to zero. In the case of Bragg diffraction, the boundary condition imposed on the 
diffracted wave is reduced to the condition that the field E^ equals zero at z = L, while 
in the case of Laue diffraction, it equals zero at 2; = 0. 

The difference between the two diffraction schemes is not merely kinematic, but radi- 
cal. Under Bragg diffraction conditions, there is a synchronous wave moving against the 
electrons of the beam, which gives rise to the internal feedback and absolute instability. In 
Laue diffraction geometry, a backward wave is absent, and as a result absolute instability 
does not evolve. It may seem that electromagnetic radiation is not generated. However, 
fluctuations, which always occur in real beams, are amplified when the beam enters the 
crystal (due to convective instability, excited in the beam). 

In analyzing multiparametric problems, to which the problem of parametric cooper- 
ative radiation refers, it is convenient to write equations (^ and ([H]) in a dimensionless 
form. This procedure enabled transferring the calculation results from one set of wave- 
lengths and beam energies to another. The analysis shows that this procedure is easy 
to perform in the case when the electrons of the beam are ultrarelativistic (the Lorentz 
factor 7^1). 

Let us note that the kinetic equation ([T]) is equivalent to the system of relativistic 
equations of motion for the momenta Pj and coordinates fj of particles, since 

j 

(the lower index j runs over all beam's electrons inside the crystal). When 7^1, these 
equations are convenient to write with the phases of charged particles (pj = k^fj — ut 
substituted for independent variables. The substitution of variables utO' — > t, uLO'/c — )■ 
L, mcujO Eq^t- / Qf, — > Eqj. {9 is the angle between the particle velocity and the wave vector 
k) then gives 

^h^2{-^-lf''ME.,e-*>), (7) 



dt'' \ dt 

dEp ^ , ^XOp , ^Xrp _ >v^ Xbj e~"^- 

dt ^^"^ dz ^ 29' "^ ^ 29'^ ~ ^ 29' N, 



dEr dEr iXOr^ , iXr J-, 



(8) 



dt ' dz 29' " 29'' 
Here the quantity Xbj = —4:TTq'nj/'muj' is determined at the moment when the jth particle 
enters the system, Uj is the corresponding electron density, and A*"; is the number of 



particles over the length /. The set of equations with boundary conditions contains four 
independent parameters: Xo,T,b/6'^, uLO"^ /c that define the geometry of the system. In 
addition to these parameters, we need to specify the beam profile. Let 



Xb = Xbo 



(1 _ ^-ct/Lf^e(ct-Li,){ct-Lb)/Lf^ 



(9) 



where the rise time Lf/c is further assumed to be equal to O.IL5/C {Lj, is the bunch length 
and O is the Heaviside function). 



3 Simulation results 



The characteristic feature of cooperative pulses is the peculiar dependence of the peak 
power on the number of particles A''^ in the bunch. When the particles are small in number, 
the radiation power monotonically increases until saturation is achieved. 
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Figure 2: Parametric radiation under Bragg diffraction conditions. Left: peak values of 
the fields |-Eo,t| plotted as a function of the bunch length. Two upper curves are plotted 
for radiation at large angels, while the two lower ones - for radiation at small angles. Solid 
curves illustrate the case when 7^ = —1; dashed curve - the case when 7,- = —0.5. Right: 
the example of a cooperative pulse [7,- = — LO]. Solid curve refers to a diffracted wave, 
dashed curve - to a direct one. 
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Figure 3: Parametric radiation under Laue diffraction conditions. Solid curve is plotted for 
radiation at large angels, while dashed curve - for radiation at small angles. [7,- = —0.5]. 
Left plot: U = 12.4. Right plot: U = 24.9. 



Let us see now how the dynamical diffraction of electromagnetic waves affects the 
cooperative radiation in crystals. Does parametric radiation possess the above described 
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Figure 4: Pulse of cooperative radiation under Laue diffraction conditions in the presence 
of noise (left) and without it (right). Solid curve shows radiation at large angles, dashed 
curve - radiation at small angles. The amplitude of white noise ratio to the average current 
density is assumed to be jrando/jav = 0.01. Computation is made for the case when the 
bunch length Lb = 49.8. 



features inherent to other types of cooperative radiation? How the Bragg diffraction case 
is different from the case of Laue diffraction? 

We shall assume that 7 = 5, Xo = Xr = 0.06, 70 = 0.98 and specify the system of 
units, where c = 1, w^^ = 1, mcud/qe = 1 and L = 24.9. 

Let us start our consideration with Bragg diffraction. The results of computation 
(Fig. 2) show that the radiation emitted at small angles reaches saturation at significantly 
smaller lengths of the bunch than that emitted at large angles. The transition from a one- 
dimensional distributed feedback (7,- = —1) to a two-dimensional one (7,- 7^—1) hardly 
causes any change in the intensity of radiation emitted at small angles to the particle 
velocity. 

Under Laue diffraction conditions, the pulse of parametric radiation has a more com- 
plicated structure due to the convective character of instability: in the absence of noise, 
the generation occurs only at the ends of the bunch of charged particles. As a result, the 
cooperative pulse possess a two-peak structure (Fig. 3). The presence of noise leads to 
an appreciable change in the pulse form: the interval between the two pulses is filled with 
a chaotic signal (Fig. 4). 



4 Conclusion 

This paper studies the features of parametric cooperative radiation emitted at both 
large and small angles to the particle velocity direction. Both Laue and Bragg cases 
are considered. It is shown that under Bragg diffraction conditions, the intensity of 
parametric radiation emitted at small angles reaches saturation at significantly smaller 
number of particles in the bunch than the intensity of radiation emitted at large angles. 

Under Laue diffraction conditions, the pulses of parametric radiation, emitted at both 
large and small angles to the particle velocity, possess a two-peak structure; each of the 
peaks is associated with the front or back end of the bunch. The presence of noise leads 
to an appreciable change in the pulse form: the interval between the two pulses is filled 
with a chaotic signal. 



A Particle-in-cell method 

The set of equations (^ and (jH]) was solved using the particle-in-cell method, which 
is widely used in plasma physics [T0| [TT] . This method implies that the solution of the 
kinetic equation is modeled using a large number of macroparticles moving along the char- 
acteristics of the kinetic equation. The current and charge densities are calculated from 
particle velocities and positions and are further used for computations of the electric field 
on a space-time mesh. The mesh values of the field are interpolated to the macroparticle 
locations; then the forces acting on macroparticles are calculated. The approach described 
here is close to the method described by I. J. Morey and C K. Birdsall in [12], which was 
used for travelling wave tube modeling. 

Let us introduce a spatial u^ = {zn = nAz,n = 0, 1, ..., rimax, '^max^^ = L} and a 
time Ut = {tg = sAt,s = 0,1,...} mesh. Specify an implicit finite-difference scheme of 
field equations ([6]) with second order accuracy in time and coordinate: 

at n = : E^+^ = 0, 

Kn - Kn ^ -^J^rn + 4J^^„+i - J^rn+2 _ ^XS)_J^s+l/2 _ ^Xr pS+1/2 

cAt ^" 2Az 2 ™ 2 °" ' 

^On ~ -^On _ 'J-'^On ~ ^-^On-l + ^On-2 *X0 ,7^+1/2 ^Xt j-,s+1/2 tS+1/2 

^At " ~^° 2Ai 2'^ ~ ^^™ ~ "^O" ' (10) 

cLTJ U *\ Ti *\ T^mnT • 



ps+l ps T^s+l/2 pS+1/2 

cAt ~ ^° 2Az 2 °" 2 ™ ^°" 



T-iS + l 7T1S ps-l-i/^ pr^-t-l/z 

^TTi ~ ^rn ^ _ -^rn+1 ~ -^rn-1 _ ]X^ ps+1/2 _ ;^ pS+1/2 

cAt ^" 2Az 2 ™ 2 °" ■ 

Let us define the source Jg„ on the right-hand side of (fTOj) . using the formula 



\^n ,.s+l/2 ^n+l ^y 



+V2 _ 2vrsinge-(^+-^V^) /^ /,^^^-^^ ^^,^^..v. 



+Eo 



s+l/2 



'11^ 



^^^^' A^ 



The contributions to each node come from the particles concentrated in the domain Zn-i < 
Zj < Zn+i- Summation in the first and second terms is made over all particles in the 
domains 2;„ < -z* < -Zn+i and z„_i < -z? < z^ respectively. The weighting factors 

s + l/2 s+1/2 

^^^ — ^j and — — ^^ are responsible for linear interpolation of the contributions to 

the node with number n that come from each particle. 

Complete the implicit difference scheme (ITO|) with the discrete analogues of the equa- 
tions of motion of macroparticles: 

P^j ~ P'^j _ 9r).R^rP^ + l/2 ifc,2j + l/2_.^^^_^^^/2)N s + 1/2 _ Pj + P j 



At 



2Q,Re(Eo7^/^e*'=^^/ -^t+dtm^^ ^ pS, 






^* •' '" . /^ + (ptr'VM,c)^ 



The field E^ at particles' locations can be found by means of linear interpolation from 
surrounding nodes: 

_ s+l/2 s+1/2 _ 

pS+1/2 _ ^n+1 - Zj j^s+1/2 ^j ~ ^n ^s+1/2 

^Oi - ^^ ^On + ^^ ^On+1 > (13) 



7 < ^'^+^/2 



Injection and extraction of particles are performed as follows: during every time step, 
we inject AA^ number of particles, whose coordinates are 



z = -Az[l + -^) , z = 0, 1, ..., AN - 1. 

Assume that each particle has a momentum p^o and a charge Q = javAt+jrandAt. Here 
jav is the average current density and jrand is a random quantity. This quantity appears 
due to shot noise of charged particles, which in high power devices can achieve quite large 
values, because the electron flow consists of ectons |13] . Since the parametric radiation is 
narrow-band by nature, and the current density spectrum jrand is a slowly varying function 
of frequency u, one can replace the real current density jrand by a random quantity 
uniformly distributed in the domain {—jrandoijrando)- The level of noise is determined 
by the fluctuation current density jrandoi^)- By analogy with the case of fluctuations of 
current, particle distribution over momenta can be introduced into the system. Particles 
are extracted from the system after they reach the coordinate Zj > L + Az. 
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